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ABSTRACT 

The existence of the Fundamental Plane imposes strong constraints on 
the structure and dynamics of elliptical galaxies, and thus contains important 
information on the processes of their formation and evolution. Here we focus 
on the relations between the Fundamental Plane thinness and tilt and the 
amount of radial orbital anisotropy: in fact, the problem of the compatibility 
between the observed thinness of the Fundamental Plane and the wide spread 
of orbital anisotropy admitted by galaxy models has been often risen. By using 
N-body simulations of galaxy models characterized by observationally moti- 
vated density profiles, and also allowing for the presence of live, massive dark 
matter halos, we explore the impact of radial orbital anisotropy and insta- 
bility on the Fundamental Plane properties. The numerical results confirm a 
previous semi-analytical finding (based on a different class of one-component 
galaxy models) : the requirement of stability matches almost exactly the thin- 
ness of the Fundamental Plane. In other words, galaxy models that are radially 
anisotropic enough to be found outside the observed Fundamental Plane (with 
their isotropic parent models lying on the Fundamental Plane) are unstable, 
and their end-products fall back on the Fundamental Plane itself. We also find 
that a systematic increase of radial orbit anisotropy with galaxy luminosity 
cannot explain by itself the whole tilt of the Fundamental Plane, becoming 
the galaxy models unstable at moderately high luminosities; at variance with 
the previous case their end-products are found well outside the Fundamen- 
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tal Plane itself. Some physical implications of these findings are discussed in 
detail. 

Key words: 
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- galaxies: kinematics and dynamics 



1 INTRODUCTION 

The Fundamental Plane (FP) of elliptical galaxies (Djorgovsky & Davies 1987; Dressier et 
al. 1987) is a scaling relation between three of their basic observational properties, namely 



the circularized effective radius {R)c = vOxA; (where Oe and he are the major and minor 
semi-axis of the effective isophotal ellipse), the central velocity dispersion ctq, and the mean 
effective surface brightness {I)c=L^/2tx{R)^ (where Lb is the luminosity of the galaxy, for 
example in the Johnson B-band). An interesting parameterization of the FP, that we adopt 
in this paper, has been introduced by Bender, Burstein & Faber (1992, hereafter BBF): 

^ l"g"g+i"g(^)-, (1) 
v2 

, _ l0gag + 2l0g(/)e-l0g(i?)e 

k2 = , (2) 

h = . (3) 

In particular, when projected on the (fci, k^) plane, the FP is seen almost edge-on and it 
is considerably thin, while the distribution of galaxies in the (/ci, /C2) plane is considerably 
broader. For example, Virgo ellipticals studied by BBF are distributed on the (fci, k^) plane 
according to the best-fit relation 

fcs = O.lSfci + 0.36 (4) 

(when adopting respectively, kpc, km s^^ and Lb© pc~^ as length, velocity and surface 
brightness units), with a very small dispersion of (y{k^) ~ 0.05 over all the range spanned by 
the data, 2.Q<ki<A.Q (and so 0.75 ^ A;3 < 1.05, see, e.g., Ciotti, Lanzoni & Renzini 1996, 
hereafter CLR96). 

By combining equations (1) and (3) with equation (4) the FP equation of BBF is then 
obtained directly in terms of the observables: the exponents are in good agreement with 
those derived (in the Johnson B-band) from a much larger galaxy sample by J0rgensen, 
Franx & Kjaergaard (1996). 
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Note that equation (4) implies that for galaxies of given luminosity Lb their effective ra- 
dius and central velocity dispersion must be strongly coupled: in fact at any fixed luminosity 
the coordinates and ki are related through definitions (1) and (3) by 

fc3 = /|A:i + ^log^, (5) 

and the slope of this relation 0.82) is different from that of the FP 0.15). As a 
consequence, in the {ki, k^) plane all galaxies with the same luminosity are located on straight 
lines significantly inclined with respect to the FP: the presence of substantial scatter in 
galactic properties from galaxy to galaxy (of similar luminosity) would destroy the thinness 
of the FP by producing a large scatter in ki and so in ^3. 

The relation between galaxy properties and the FP can be expressed in a quantitative 
way under the reasonable assumption that present day ellipticals are virialized systems. We 
write the virial theorem as 

where T,, = M^/L-q is the galaxy stellar mass-to-light ratio (for example in Blue solar 
units), and is a dimensionless factor depending on the stellar density profile, internal 
dynamics, dark matter amount and distribution, and, for non-spherical galaxies, on their 
relative orientation with respect to the observer's line-of-sight (see, e.g., Ciotti 1997); in 
addition, K\ depends also on the observing aperture adopted to derive ctq. 

Equations (4) and (6) imply that, for galaxies belonging to the FP, the quantity T^/Ky 
is a very well defined function of two0 of the three observables Lb, ctq, and {R)c- In the par- 
ticular case of adopting the numerical coefficients of equation (4), measuring Lb in IO^^Lb©, 
and taking into account that ki = \og{GT ^L^/ Ky) / \^ and k^ = log(27rG'T^,/_K"v)/-\/3, 

^ ^ 1.12 X Lb°-^^ (7) 

Ay 

where the quantity T^,/i^v is characterized by a scatter of ~ 20% due to its relation with 
^3 (however, additional considerations reduce this figure to ~ 12%, see Renzini & Ciotti 
1993). As a consequence, any departure from the relation dictated by equation (7) will move 
a galaxy away from the FP. We recall here that the dependence of k^ on ki, as given by 

* In principle one could find virialized galaxies everywhere in the three-dimensional observational (Lb,(to, {R)c) space. From 
this point of view the existence of the FP is related to the virial theorem as the HR diagram is related to hydrostatic equilibrium: 
the useful information that one derives is not about the equilibrium equations, but the physics of the objects involved. 
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equation (4) and responsible for the luminosity dependence of the ratio T^-jKy, is commonly 
known as the FP "tilt" . 

The simple analysis presented above shows that the two properties of thinness and tilt 
of the FP are deeply connected with the present-day structure and dynamics of ellipticals 
(hereafter, Es), and, as a consequence, with their formation and evolution history: the very 
existence of the FP [as well as of the other tight scaling relations revealing the remarkable 
homogeneity of Es, such as the Mbh — o"o (Gebhardt et al. 2000; Ferrarese & Merritt 2000), 
the Mgg-CTo (Bender, Burstein & Faber, 1993 and references therein) and the color-magnitude 
(Bower, Lucey & Ellis, 1992) relations] imposes strong constraints on the different formation 
and evolutionary scenarios proposed for Es (i.e., dissipationless merging, monolithic collapse, 
or a combination of the two; see, e.g., Ciotti & van Albada 2001). 

Among the various galaxy properties in principle able to destroy the FP thinness (as a 
consequence of a substantial variation at fixed galaxy luminosity), one of the most "effective" 
is certainly orbital anisotropy (de Zeeuw & Franx 1991). In fact, galaxy models are com- 
monly believed to be able to sustain a large spread of orbital anisotropies and it is also well 
known that radial orbital anisotropy can produce very high central velocity dispersion values, 
and correspondingly low values of Ky, thus substantially violating equation (7). A natural 
question to be addressed is then what physical principle or evolutionary process limits the 
range of orbital anisotropies shown by real galaxies. Ciotti & Lanzoni (1997, hereafter CL97), 
using one-component, radially anisotropic Sersic (1968) models, and a semi-analytical in- 
vestigation based on the Fridman & Polyachenko (1984, hereafter FP84) stability indicator, 
suggested the possibility that radial orbit instability could be the limiting factor of the FP 
thickness. In practice, CL97 found indications that galaxy models sufficiently anisotropic to 
be outside the FP observed thickness (when their parent isotropic model was assumed to lie 
on the FP) were unstable. Clearly, this preliminary indication requires a confirmation with 
the aid of numerical simulations and more realistic galaxy models (for example, allowing 
for the presence of live dark matter halos). Also, a question naturally associated with that 
above is the determination of the position of the end-products of unstable initial conditions 
in the space of the observables. One of the aims of this work is indeed to answer these two 
questions by numerical simulations of one and two-component radially anisotropic galaxy 
models. 

Orbital anisotropy is not only related to the problem of the FP thinness but it is also 
one of the candidates that have been proposed to explain the origin of the FP tilt (CLR96; 
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CL97). If this were the case, the amount of radial anisotropy in the velocity dispersion 
tensor should increase with galaxy luminosity, as can be seen from equation (7) under the 
assumption of a constant T* and of structural homology]^ over the whole FP plane. In 
other words, in this scenario the FP tilt would be produced by a dynamical non-homology 
due to anisotropy (note that dynamical non-homology may well coexist with structural 
homology, but the converse is in general not true). Many interesting questions are raised 
by the scenario depicted above: for example, under the assumption that an isotropic galaxy 
of given luminosity lies on the FP, how far can the derived structurally homologous but 
radially anisotropic models climb over the FP before the onset of radial orbit instability? 
In addition, what happens to the end-products of the unstable models? Will they remain 
near the FP? In this paper we try to address also these questions with the aid of N-body 
numerical simulations. 

Strictly related to the clarification of the interplay between orbital anisotropy and the 
FP tilt, is the possibility to obtain some clues on the formation processes of Es. In fact it is 
trivial to prove that, as consequence of the virial theorem and the conservation of the total 
energy, in the merging of two galaxies with masses Mi and M2 and virial velocity dispersions 
o"v,i and o"v,2? the virial velocity dispersion of the resulting galaxy is given by 
2 _ Miali + M2al^, 

(by definition = 2T/M, where T is the total kinetic energy of the galaxy). For simplicity 
in the formula above we considered, in the initial conditions, a negligible energy of the galaxy 
pair when compared to the other energies involved in the process, and no significant mass loss 
from the resulting system. From equation (8) it follows that 0-^,1+2 < max(crv 1, o"v,2), i-e., the 
virial velocity dispersion cannot increase in a merging process of the kind described above. 
On the other hand, the FP (or the less tight Faber-Jackson relation; Faber & Jackson 1976) 
indicates that the projected, central velocity dispersion increases with galaxy luminosity. 
Then, in the dissipationless merging scenario the FP tilt can be produced only by structural 
and/or dynamical non-homology, since the relation between central and virial velocity dis- 
persion depends on the structure and dynamics of the system: in particular we will discuss 
here the second possibility, with regard to an increase in radial orbital anisotropy with galaxy 
luminosity. Note that an increase in the radial orbit amount has been claimed in the past as 

t With "structural homology" we mean that all the structural galaxy properties (e.g., the stellar and dark halo density profiles, 
the ratio of their scale— lengths and masses, and so on) do not depend on Lb- 
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a natural by-product of galaxy merging, and also some observations have been interpreted 
in this way (Bender 1988; see also Naab, Burkert & Hernquist 1999, and references therein). 
With the aid of the explored numerical models we will try to obtain some qualitative insight 
in this problem: it is however clear that the results should be considered at the best qual- 
itative indications, and that a firm answer about the role of merging in producing the FP 
tilt can be obtained only with N-body numerical simulations of merging galaxies (see, e.g., 
Capelato, de Carvalho & Carlberg 1995). 

Summarizing, the aims of this work are the following. For what concerns the FP thickness 
problem, we investigate the role of radial orbit instability as a factor regulating the amount 
of radial anisotropy for galaxies of given luminosity, and the position, relative to the FP, of 
the end-products of radially unstable anisotropic models. In addition, we determine whether 
it is possible to reproduce the whole FP tilt with a systematic variation of radial anisotropy 
with luminosity (using both stable and unstable initial conditions), and what is the fate 
of unstable models initially forced to lie on the FP. This paper is organized as follows. 
In Section 2 we describe the basic structural and dynamical properties of the investigated 
models; a short description of the codes used for the numerical simulations is also given. 
In Section 3 we present the results and their impact on the FP thinness problem, and in 
Section 4, the results are discussed focusing on the origin of the FP tilt. Finally, in Section 
5 the main conclusions are summarized. 

2 MODELS 

2.1 Initial conditions 

2.1.1 Structural and dynamical properties 

As initial conditions for the N-body simulations we use spherically symmetric one-component 
7-models (Dehnen 1993; Tremaine et al. 1994) and two-component (71,72) models (Ciotti 
1996, 1999). The density, mass, and (relative) potential profiles of the stellar component are 
given by 



3-7 M,re 
An r-^ir-c + ry-^ 



(0 < 7 < 3), 



(9) 




(10) 




(11) 



© 2001 RAS, MNRAS 000, 



Radial orbital anisotropy and the Fundamental Plane of elliptical galaxies 7 



(7 = 2), 



(12) 



where is the total stellar mass; the 7=1 and 7 = 2 cases correspond to Hernquist (1990) 
and Jaffe (1983) density distributions, two reasonable approximations (when projected) of 
the i?^/^ law (de Vaucouleur 1948). In the two-component models, the dark matter (DM) 
halo is described by ph, and profiles of the same family of equations (9)-(12), where 
now r\i = f3rc and Mh = pM^,. 

Radial anisotropy in the stellar orbital distribution is introduced by using the Osipkov- 
Merritt (OM) parameterization (Osipkov 1979; Merritt 1985). The distribution function 
(DF) of the stellar component is then given by 



The variable Q is defined a.^ Q = £ — L^/2r^, where the relative (positive) energy is given 
by £^ = \E't — v"^ I'^i v is the modulus of the velocity vector, the relative total potential 
is \E't = + \I'h, L is the angular momentum modulus per unit mass, and f*{Q) = for 
Q < 0. The quantity ra is the so-called "anisotropy radius" : for r ^ the velocity dispersion 
tensor is mainly radially anisotropic, while for r <^ ra the tensor is nearly isotropic. Isotropy 
is realized at the model center, independently of the value of r^; in the limit Ta 00, Q = S 
and the velocity dispersion tensor becomes globally isotropic. 

In order to reduce the dimensionality of the parameter space, the orbital distribution of 
the DM halos is assumed isotropic in all our simulations; as a consequence the DF for the 
DM halo component, fh{£), is given by equation (13) where Q = £ and Ph(^) substitutes 

2.1.2 Physical scales 

According to the given definitions, from the structural and dynamical point of view the one- 
component models are completely determined by four quantities: the two physical scales 
and Tc, and the two dimensionless parameters 7 and Sa = Ta/rc. 

In the numerical simulations M^, Tc and Tdyn are adopted as mass, length and time scales. 
Tdyn is the half-mass dynamical time, defined as 




(13) 



where 




(14) 
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where pu = 3M=K/87rrfj is the mean density inside the half-mass radius tm and 

F(7) = (2^-l)"^ (16) 

Note that at fixed M^, and Vc the half-mass dynamical time depends strongly on 7. For 
example, F(0)/F(1) ~ 2.0 and F(l)/F(2) ~ 3.8. Finally, the velocity scale is given by 



In the case of two-component galaxies we limit our present study to (1,1) models, i.e., 
to two-component Hernquist models (Ciotti 1996). In this way the DM halo is similar to 
Navarro, Frenk & White (1996) profiles. (1,1) models are characterized by five quantities, 
the two physical scales M^, and r^, and the three dimensionless parameters Sa, p and (3. Their 
Tdyn depends on /i and /5: from definition (15), where now pu = 3(1 + p)M^/ (87rr|[) and Tm 
is the half-mass radius of the total (stellar plus dark) density distribution, it results that 



where F = x^/^(/i, /5)/vT+7'') and x is the solution of 



Finally note that, once the dimensionless numbers Sa and 7 (for one-component mod- 
els) and Sa, p, and (3 (for two-component models) are fixed, the results of the numerical 
simulations can be rescaled for arbitrary values of M^,, r^, and Lb (or T^,). In particular, in 
order to transform the results of the numerical simulations in physical units we express 
in 1O^°M0, rc in kpc and Lb in IO^^Lb©. 

2.1.3 Numerical realization of the initial conditions 

In order to arrange the initial conditions for the numerical simulations, we distribute N 
particles (A^ = 32768 and = 131072 for one and two-component models, respectively) by 
using spherical coordinates (r, t?, for positions, and Vt-,v^,v^, for velocities). The radial 
coordinate is assigned by inverting equation (10) and choosing a sampling suited to resolve 
the core of the distribution, while angular coordinates are given randomly. The velocity 
vector of each particle is assigned by using the von Neumann rejection method. This method 
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can be easily applied to OM systems by introducing the dimensionless vector of components 
{ui,U2,U3), related to {vr,v^,V:^) by 



Vr = ^JZWtUi, = U2, = Ms, (20) 

where 



1 + 'i (21) 

a 



from this choice < m < 1, where u = ywf + U2 + u^. The vector (ui,U2,M3) is assigned 
to each particle by using the rejection method, after computing the value of the DF at 
the required Q = {1 — n^)^T, by numerically evaluating equation (13). The components 
of the velocity vector {v-r:,v^,v^) are then recovered according to equation (20). The initial 
conditions so obtained are then compared to the expected analytical density distribution, 
finding, as a rule, an agreement of better than 2% over the radial interval containing 0.99 of 
the total mass of the theoretical model. 

2.2 The numerical codes 

For our simulations we used the Barnes & Hut (1986) TREECODE (in a version made pub- 
licly available by Hernquist, 1987) and the Springel, Yoshida & White (2000) parallel version 
of GADGET (adapted to run on the 32 processors Cray T3E at CINECA). In particular, 
the one-component simulations were run on a Alpha workstation by using the TREECODE 
with = 32768; for some of them we used also A^ = 65536. The quadrupole correction in the 
cell-particle force calculation was always applied. For the two-component models we used 
instead GADGET with A^ = 131072; in addition, some one-component models were also 
tested with A^ = 131072 using the same code. We found very good agreement between the 
test simulations performed with both the codes, with the basic properties of the simulated 
systems nearly independent of the adopted number of particles. 

From the numerical point of view, TREECODE simulations are characterized by three 
parameters, namely the opening angle 6, the softening parameter e (i.e., the softening length 
expressed in units of Tc) and the time step At. Following Hernquist (1987) and Barnes & 
Hut (1989), we assume 6 = 0.8, a good compromise between conservation of total energy 
and computational time. The choice of At and e requires some care; in fact these two 
parameters are strongly coupled and, in order to maintain the same accuracy in the force 
evaluation. At must be reduced if e is reduced (Barnes & Hut 1989). In addition, the 
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"optimal" value of the softening length is strongly dependent on and on the specific density 
distribution profile (Merritt 1996; Athanassoula et al. 2000; Dehnen 2001). To choose the 
values of these parameters, we performed some simulations of isotropic 7-models over 100 
Tdyn, checking the total energy E and virial ratio {V = 2T/W) conservation. On the basis 
of these tests we adopt At = Tdyn/100, and e = (0.072, 0.030, 0.016) for 7 = (0, 1, 2) models, 
respectively. With this choice we obtain \AE/E\ < 1.5% and \AV/V\ < 2%. GADGET 
simulations depend on five parameters: the cell-opening parameter a, the minimum and 
the maximum time step Atmin and Atmax, the time-step tolerance parameter atoi, and the 
softening parameter e (Springel et al. 2000). Fiducial values of the parameters were fixed 
after running a few one-component (7 = 1) test simulations. In particular, with a = 0.02, 
Atinin = 0, Atmax = Tdyn/lOO, atoi = 0.05, and e = 0.016 we obtain \AE/E\ < 0.5% and 
|Ay/y| < 1%. The value of the first four parameters were used also in the two-component 
simulations. The softening length e was instead determined case by case, by considering the 
core radius length of the more concentrated component. 

In order to determine the "observational" properties of the end-products of the numerical 
simulations, and place them in the (fci, k2, k^) space, we measured their effective radius (i?)e, 
central velocity dispersion ctq and mean effective surface brightness (J)e for several projection 
angles. The technical details of the adopted procedures are given in the Appendix, together 
with a brief discussion on the discreteness effects on the derived values of the "observational" 
properties of the models. A point of interest is the choice of the adopted "aperture radius" 
for the measure of ctq, fixed in this work to {R)c/8 in order to match the observational 
procedure at low redshifts (see, e.g., J0rgensen et al. 1996). Indeed, it is well known that 
isotropic 7-models may present a projected velocity dispersion profile decreasing toward the 
center (in contrast to what happens in the vast majority of elliptical galaxies): for example, 
the projected velocity dispersion of Hernquist models peaks approximately at Re/5, while 
the profile for (isotropic) Jaffe models is monotonically decreasing. As a consequence, the 
position of the initial conditions and of the end-products in the k space could depend on the 
adopted aperture radius used to measure (Jq. In order to assess the effect of this choice on 
our conclusions, we analyzed the results of the simulations also by using an aperture radius 
of (i?)e/4, and we found very good agreement with the results obtained with the aperture 
{R)e/8: therefore, we present here only these last results. Note however that, by virtue of 
the projected virial theorem (see, e.g., Ciotti 1994), a totally different scenario from that 
explored in this paper would arise in the limiting case of (Jq measured over the whole galaxy: 
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Table 1. The dimensionless coefficient i<'v(7, Sa) for one-component models, as obtained from equation (6). See Section 3 for 
the definitions of Sac and Sass- 

7 Kv{-y,oo) Ky{'y,Sas) -R'v(7iSac) 

6.7 5.9 2.8 

1 5.9 5.3 2.7 

2 4.8 4.4 2.3 

anisotropy would play no role at all in determining the position of initial conditions in the k 
space, and any difference between initial and final o"o would be only due to projection effects 
associated with loss of spherical symmetry. These aperture effects could be important when 
studing observationally the FP at intermediate redshift (see, e.g., van Dokkum & Franx 
1996, Bender et al. 1998, Pahre, Djorgovski & de Carvalho 1998, Treu et al. 1999), or locally 
by using large apertures (see, e.g., Graham & CoUess 1997). 

3 ORBITAL ANISOTROPY AND THE FP THINNESS 

As already pointed out in the Introduction, the 1-sigma dispersion of the observational 
data around the best-fit relation (4) is surprisingly small and nearly constant over all the 
observed range in ki, with a{ks) ~ 0.05. We investigate here the constraints imposed by this 
tightness on the amount of radial orbital anisotropy for the set of one and two-component 
galaxy models described in Section 2. We start by fixing the values for the dimensionless 
parameters 7 (for one-component models) and fi and (3 (for two-component models); we 
also assume global isotropy, i.e., Sa 00, and in this way the quantity Ky is uniquely 
determined. These globally isotropic models (that we call parent models) are then placed on 
the FP by assigning of the pair (T^,, Lb) so that equation (7) is verified. From each of these 
parent models lying on the FP we then generate a family of CM radially anisotropic models 
by decreasing s^, while maintaining fixed all the other model parameters. Correspondingly 
Ky decreases (as can be seen from Table 1, in the case of one-component models and 
for representative values of the parameter Sa), ki increases, and so does k^, according to 
equation (5): for sufficiently small values of Sa the members of each family are found outside 
the observed thickness of the FP. Note that by an appropriate choice of the pair (T^,,Lb) 
each parent isotropic model can be placed at arbitrary positions over the best-fit line (4), 
and so the results of the numerical simulations (after a rescaling to T^, and Lb) are the same 
everywhere on the FP. In addition, since alk^) is constant over the whole observational range 
spanned by ki, the conclusions obtained from each family of models are also independent of 
the position of the parent galaxy on the FP. 
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Obviously all isotropic models discussed in this paper are stable (see, e.g., Binney & 
Tremaine 1987; Ciotti 1996), while for each family of models a critical value Sas for stability 
exists such that the initial conditions characterized by Sa < s^s describe radially unstable 
configurations. From the point of view of the present discussion, the critical value Sas cor- 
responds, through irv('5as); to the maximum distance that a stable model can have from 
the FP, where its parent isotropic model lies at (say) {kf",k^^°). Clearly, initial conditions 
describing unstable models can be placed at larger distances from the FP: these initial con- 
ditions will evolve with time, and their representative points in the space of observables will 
also evolve with time, up to virialization. In particular, the coordinates ki{t) and k^it) will 
evolve with time moving on the line described by equation (5). Note that the maximum 
distance from the FP at which unstable models can be placed is in general finite: in fact, 
the anisotropy radius of all physically acceptable (stable and unstable) galaxy models must 
satisfy the inequality Sa > Sac, where Sac < Sas is the (dimensionless) critical anisotropy 
radius for consistency (i.e., the anisotropy limit for initial states with a nowhere negative 
DF); for a study of s^c in 7-models and in (71,72) models see CaroUo, de Zeeuw & van der 
Marel (1995) and Ciotti (1996, 1999). 

FP84 argued that a quantitative indication on the maximum amount of radial orbits 
sustainable by a specific density profile is given by the stability parameter ^ = 2T^/Tt, 
where Tr and Tt = + are the radial and tangential component of the kinetic energy 
tensor, respectively. From its definition ^ 1 for Sa ^ 00 (globally isotropic models), while 
^ — i> 00 for Sa — * (fully radially anisotropic models). The fiducial value indicated by FP84 
as a boundary between stable and radially unstable systems is — 1-7. Unfortunately, the 
reliability of such indicator is not well understood, and indications exist of a significant de- 
pendence on the particular density profile under scrutiny (see, e.g., Merritt & Aguilar 1985; 
Bertin & Stiavelli 1989; Saha 1991, 1992; Bertin et al. 1994; Meza & Zamorano 1997). In 
any case, CL97 used this value to determine Sas for one-component Sersic (1968) models by 
solving the associated Jeans equations, and from this value they determined the maximum 
distance of stable models from the FP: all models characterized by ^ < where found inside 
the observed thickness of the FP. This finding can be considered at the best a qualitative 
indication, considering the uncertainties associated to the exact value of .^s and to its de- 
pendence on the specific density profile adopted, and the need of numerical simulations is 
clear. 

With the aid of N-body simulations in the following two Sections we investigate how 
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distant from the FP stable models of various families can be placed, by increasing their radial 
orbital anisotropy. The logically related question of what happens to the end-products of 
the unstable (but physically consistent) initial conditions is also addressed. 

Due to its central role in the following discussion, it is important to quantify the concept 
of "distance" of a galaxy model from the FP. In general, we define distance of a point (fci, k^) 
from the FP the quantity 6ks = j/cs— O.lS/ci— 0.36|, i.e., the distance at fixed ki from the point 
and the FP itself: this quantity naturally compares with cr(/c3). Unfortunately, as already 
discussed in the Introduction, galaxy models of fixed luminosity move along inclined lines in 
the {ki, ks) space and so in the present exploration \k3 — kf°\ is not the distance from the FP. 
The relation between 6k^ and | ^3 — is however of immediate determination: in fact, from 
equation (5) k^ = kf° + ^J2/3{kl — kf°), and from the assumption that the parent isotropic 
models are placed on the FP one obtains Sk^ = (1 — 0.15-^3/2)| /cs — ~ 0.816|A;3 — kf°\. 

3.1 One— component models 

In order to answer to the questions outlined above, we performed a set of 21 simulations 
of one-component 7-models with 7 = (0, 1, 2), evaluating numerically for each of the three 
families the critical value for stability Sas- We recall here that in our simulations the onset 
of instability is just due to numerical noise produced by discreteness effects in the initial 
conditions, and that, at most, the simulations are interrupted after 100 T^yn- As a rule, in 
order to determine if a given model is unstable, we found useful to check its departures from 
spherical symmetry by monitoring the evolution of its intrinsic axis ratios c/ a and b/ a (where 
a, b and c are the longest, intermediate and shortest axis of its inertia ellipsoid associated 
to "bona fide" bound particles, see Appendix): we found that numerical uncertainties (due 
to the finite number of particles) of these ratios never exceed 5%. According to this choice 
we define unstable the models for which the minimum c/a over 100 T^yn is smaller than a 
fiducial threshold value, 0.95 (the horizontal dashed line in Fig. 1). Also this value has been 
obtained by analyzing the fluctuations - due to the finite number of particles - of c/a shown 
by the numerical realizations of the isotropic (stable) parent 7-models. As expected, we 
found that an exact determination of s^s for a given density profile is not straightforward: in 
fact, while for strongly anisotropic initial conditions the onset of instability is apparent and 
the numerical models settle down into a final equilibrium configuration in a few dynamical 
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Figure 1. Left: final axis ratio (c/a)^^ vs. the stability parameter ^ = 2Tr/Tt of initial conditions, for all the computed models. 
One— component models are represented by full symbols (7 = 0, circles; 7 = 1, squares; 7 = 2, triangles). Two— component (1,1) 
models are represented by empty symbols: squares correspond to concentrated halos (/3 = 0.5), circles to diffuse halos (/3 = 2). 
Crosses indicate models with massive halos (fi = 3). Right: final axis ratio {c/a)^^ vs. r^/rM for the same models shown in left 
panel. The vertical dashed line is located at = 1-7 (FP84). Models above the horizontal dashed line are stable (see the text 
for a discussion). 

times, for nearly stable initial conditions the instability can be characterized by very slow 
growth rates and its effects become evident even after 30 Tdyn (see also Bertin et al. 1994). 

The result of the simulations are summarized in Fig. 1, where we plot for all models the 
final value of the axis ratio, (c/a)fin, as a function of ^ (left panel) and r^/rM (right panel) of 
initial conditions (we recall here that tm is the spatial half-mass radius). A first result is that 
for one-component 7-models (full symbols) the ^ critical value is in the range 1.6 ^.^s ~ 1-8: 
for ^ < 1.6 all models were found stable up to 100 Tdyn, while for ^ > 1.8 all models present 
clear evolution on time-scales shorter than 30 T^yn- This range for is compatible with the 
value 1.7, reported by FP84 and used in CL97, and with the results of Bertin et al. (1994), 
who estimated .^s — 1-58 for the family of "/oo models"; Meza & Zamorano (1997) found 
instead a higher threshold value for stability (^s ~ 2.3). When expressed in terms of ra/rM 
the critical value for stability is found in the range 0.6 ^ (ra/rM)s ~ 0.8. Finally, the stability 



© 2001 RAS, MNRAS 000, 



Radial orbital anisotropy and the Fundamental Plane of elliptical galaxies 15 



limits expressed in terms of the (dimensionless) critical anisotropy radius (the quantity 
of direct interest in this work) are given by Sas = (3.0,1.8,0.6) for 7 = (0,1,2) models, 
respectively; for comparison the critical anisotropy radius for consistency is Sac= (0.5,0.2,0.0; 
Ciotti 1999). 

For what concerns the internal structure of the end-products of unstable initial conditions 
we found (in accordance with previous results, see, e.g., Merritt & Aguilar 1985, Stiavelli 
& Sparke 1991) that they are in general prolate systems, with axis ratios in the range 
0.3 ^ (c/a)fin ^ 1, consistent with the ellipticities of the observed galaxies. Only the most 
anisotropic models, near the consistency limit (sa — Sac), form a triaxial bar. From Fig. 1 it 
is apparent that (c/a)fin is strongly anti-correlated with ^ (and so correlated with ri^/ru), 
and this in a way essentially independent of the value of 7; a similar decrease of (c/a)fin with 
Sa was also found in the numerical simulations of Meza & Zamorano (1997). 

In order to better compare the end-products of unstable models with real galaxies, we 
also fitted their projected mass density profiles with the widely used Sersic (1968) i?^/™ law: 



where b{m) ~ 2m — 1/3 + 4/405m (Ciotti & Bertin 1999). In Fig. 2 we plot the Sersic best 
fit parameter m as a function of ^ for a small, but representative, set of models. As for real 
galaxies, we found a significant dependence of m on the adopted radial range over which 
the fit is performed, while the value of m is not very sensitive to the specific fitting method 
adopted (see, e.g., Bertin, Ciotti & Del Principe 2001). For example, in the radial range 
0.1 ^ -R/(-R)e ^ 4 we found l^m<5, with average residuals between the data and the fits 
(A/i) ~ 0.02 0.14magarcsec~^. Clearly, the fitted quantities m and {R)e depend on the 
relative orientation of the line-of-sight and of the end-products of the simulations: with the 
vertical bars in Fig. 2 we indicate the range of values spanned by m when projecting the 
final states along the short and long axis of their inertia ellipsoids. As an example of fit for 
a specific model, in Fig. 3 we show the data relative to the initial conditions of an unstable 
7 = 1 model, and to the two projections of its end-products. 

Having determined for each family of galaxy models the critical anisotropy radius, we 
can now proceed to check how distant from the FP can stable models be placed, and where 
the end-products of unstable initial conditions are found. Of course, in this second case the 
coordinates in fc-space depend on the line-of-sight orientation with respect to the density 
distribution of the end-products, and the dependence is expected to be stronger for smaller 
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Figure 2. Sersic best fit parameter m vs. stability parameter ^ for a subset of the computed models. Symbols are the same 
as in Fig. 1. The three points at ^ = 1 give the value of m for the isotropic (stable and spherically symmetric) parent galaxies, 
while vertical bars show the maximum and minimum m for each model due to the relative orientation of the end-product and 
the line-of-sight. 

values of {c/a)^^^: to any unstable initial condition corresponds, in the /c-space, a set of 
points. Due to the fact pointed out above that the properties of the models here investigated 
do not depend on the specific position of the parent galaxy on the FP, in Fig. 4 we plot 
the obtained results in a coordinate system that reflects this property, and allows for an 
immediate visualization of the most important consequences derived from the simulations; 
in Fig. 5 the behavior of a representative set of models is shown in the standard {ki, k2, k^) 
space. In the horizontal axis of Fig. 4 we plot the displacement of the anisotropic initial 
conditions, measured by their fcg, with respect to the isotropic parent galaxy, placed at 
/c3°: as a consequence, initial conditions with larger \k\ — k^^°\ are characterized by larger 
amounts of radial orbital anisotropy. As already discussed in the Introduction, these initial 
conditions are placed in the {ki, k^) space along the lines given by equation (5). For example, 
in the upper panel of Fig. 5 the initial conditions Ai, Bi, and Ci correspond to the isotropic 
parent galaxies placed at points A, _B, and C (by assuming T^, = 5 and determining L-q 
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Figure 3. Circularized surface briglitness profiles (vertical bars) of the end— product of the unstable one— component 7 = 1 
model with Sa, = 0.7. Re is the effective radius of the (spherically symmetric) initial condition (solid dots), while {R)c = V Qcbc 
is the circularized effective radius. Solid lines represent the best fit Sersic models of the end-product projections along its inertia 
ellipsoid minor (c) and major (a) axis. 

from equation [7]), while in the lower panel the same parent galaxies and initial conditions 
are shown in the {ki, /C2) space, connected by dotted lines. 

On the vertical axis of Fig. 4 we plot the quantity — kf" corresponding to the end- 
product of each explored initial condition: if k^ = kf" then the model has "fallen back" on 
the FP. We also recall here that Sk^ ~ 0.816|/c3 — k^^°\: in Fig. 4 the two horizontal dashed 
lines correspond to the FP thickness cr(A;3)/0.816 ^ 0.0613, and so points inside this strip 
represent models consistent with the observed thickness of the FP. Finally, as an obvious 
consequence of the choice of the coordinate axes in Fig. 4, note that all the initial conditions 
(as for example models Ai, Bi, and Ci) are located at t = on the dotted line (with slope 
equal to 1). This means that this line is also the locus of stable models, while all the parent 
galaxies (as for example models A, B, and C of Fig. 5) lie at the point (0, 0). A few stable 
anisotropic models may in fact be seen, as single solid points, on this line at low anisotropy 
values. 
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Figure 4. Final vs. initial kg for all galaxy models, measured with respect to fejf°, the fes coordinate of their isotropic parent 
galaxy. The horizontal dashed lines correspond to S^ks) = (7(^3), while the dotted line fca = fcg is the locus of the initial 
conditions and of the stable models. Symbols are the same as in Fig. 1. See the text for the explanation. 



By increasing the amount of radial anisotropy the initial conditions move along the dotted 
line, and when they reach the critical value of the anisotropy radius they become unstable 
and rearrange their density profile and internal dynamics in a new, stable configuration. 
As we have seen, these end-products are strongly asymmetric, and so their representative 
points span a range of values as a function on the line-of-sight orientation. In Fig. 4 the 
vertical lines show the importance of this projection effect: note that, in general, the length 
of these segments is considerably smaller than the FP thickness. 

The first result that can be obtained by inspection of Fig. 4 is the fact that, for one- 
component 7-models (solid symbols), radial orbit instability becomes effective for initial 
conditions inside the FP thickness, thus providing strong support to what found by CL97 
for i?^/"^ models. On the contrary, by considering as anisotropy limitation the very basic 
requirement of model consistency only, it is apparent from Fig. 4 that physically acceptable 
initial conditions could be placed at a distance from the FP substantially larger than ^(/cs). 
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up to Sk^ ~ 2^3cr(A;3), the well known problem motivating this work. Models Ai, Bi, and Ci 
in Fig. 5 are just three examples of such models. The second result is that all one-component, 
radially unstable models fall back on the FP: in other words, not only the FP thickness is 
nicely related to stability, hut the FP itself acts as an "attractor" for the end-products of 
radially unstable systems when their parent galaxies lie on it. 

The same results can be illustrated in a more direct way by using Fig. 5, where we plot 
the positions of the end-products of unstable models in the k space. As described above, to 
each final configuration derived from unstable initial conditions corresponds a set of points, 
depending on their relative orientation with respect to the observer's line-of-sight. Due 
to the fact that the total luminosity of each model is obviously conserved by projection, 
these sets in the {ki, k^) space are actually segments of the straight line given by equation 
(5), as shown by Fig. 5 where the end-products of the initial conditions Ai, and Ci 
can be immediately recognized. On the contrary, in the (/ci,fc2) plane the end-products 
are distributed, as a consequence of projection along different angles, on two-dimensional 
regions. This is due to the fact that no 1-1 relation similar to equation (5) exists between ki 
and /c2, because also the value of {R)e enters explicitly [k2 = log{T ^ L'^ / Ky{R)^) + const]. 
For this reason for each end-product we plot several positions obtained with random viewing 
angles. Note however that the displacement in the {ki, ^2) space occurs mainly along the ^2 
coordinate, and this is due to the steepening of the profile (Fig. 2), i.e., to the increase of 
(/)e. In any case, the models remain well inside the populated zone of the {ki, /C2) plane. 

3.2 Two— component models 

We now present the results of the numerical simulations (8) of two-component (1,1) galaxy 
models. In fact, although indications exist that the onset of radial orbital anisotropy is 
not strongly affected by the presence of a massive DM halo (see, e.g., Stiavelli & Sparke 
1991, Ciotti 1996), in principle its presence could significantly modify the structural and 
dynamical properties of the end-products of unstable initial conditions, and so alter the 
findings obtained for one-component models described in Section 3.1. 

Unfortunately, due to the dimension of the parameter space, we are not in the position 
to determine, as for one-component models, even a fiducial threshold for stability, and so 
we limit our study to the behavior of some representative models. In particular, we consider 
the following cases: ;U = 1 ("light" halo), /i = 3 ("massive" halo), P = 0.5 ("concentrated" 
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Figure 5. Top: fcs vs. fci for a representative set of the one and two— component models shown in Fig. 3 and Fig. 4: symbols 
are the same. The solid line represents the FP relation (equation [4]) with its observed dispersion (dashed lines). The stellar 
mass-to— light ratio is fixed to = 5 and = 2.7, for one and two-component models, respectively. Bottom: k2 vs. fci for the 
same models plotted in the upper panel. The dashed lines define the region, as given by BBF, where real galaxies are found. 
The arrow shows the direction followed by initial conditions with increasing radial anisotropy. 
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halo) and (3 = 2 ("diffuse" halo), where for each of the four possible combinations we fix the 
anisotropy radius to Sa = 0.3 and Sa = 0.7, two values corresponding to strongly unstable 
one-component 7 = 1 models. In all the simulations we use live DM halos, and in order 
to have equal mass particles for "stars" and "dark matter" we adopt A^h = 98304 and 
= 32768 for /i = 3, and = = 65536 for = 1. 

In analogy with the one-component case, the quantities Tr and Tt, entering the definition 
of ^ = 2Tr/Tt, are now the radial and the tangential component of the total (stellar and halo) 
kinetic energy tensor, respectively. This choice seems the natural one in the case of a live 
DM halo, when from a dynamical point of view the galaxy should be considered as a whole; 
but certainly other choices (for example, by using in the ^ definition the kinetic energies 
of the stellar component only) could be equally well motivated. Clearly, the determination 
of the observational properties of the end-products is based on the analysis of their stellar 
component only. 

As can be seen from Fig. 1, where the ellipticity of the stellar component of (1,1) models 
is represented by empty symbols, the basic trend of (c/a)fin with ^ is similar to that of one- 
component models: the end-products are mostly prolate systems, with axis ratio (c/a)fin in 
the same range of that of the one-component models. In general, however, when considering 
one and two-component initial conditions with the same ^, the final stellar distribution 
remains more spherical in the cases with DM than in the one-component cases. Moreover, 
models with massive DM halos (empty symbols with crosses in Figs. 1, 2, 4 and 5) remain 
more spherical than the corresponding models with light halos and (approximately) the 
same ^. This means that for the explored two-component models, at variance with the 
one-component cases, the parameter ^ is not well correlated with the final axis ratio. As 
can be seen from Fig. 2, the best fit parameter m of the projected stellar distribution of 
the end-products of (1,1) unstable models remains limited to values m<3, the same range 
covered by one-component 7 = 1 models. We also found that the final shape of the DM 
halos remains nearly spherical, with 0.88 ^ (c/a)fin ^ 1. 

The observational properties of the end-products of the two-component models are 
illustrated in Figs. 4 and 5 where a comparison with the one-component models can be 
easily made. In particular, as in the one-component case, we found that the FP thickness still 
nicely separates stable from unstable models, independently of the amount and distribution 
of DM. However, the final position in the k space of the end-products depends on the amount 
of DM: as expected two-component models with a light halo (quite independently of its 
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concentration) are very similar to one-component models, while models with massive DM 
halos are quite different. In particular from Fig. 4 it is apparent that massive halos prevent 
the models from falling back on the FP and this effect is stronger for more concentrated 
halos. This is shown in Fig. 5 (upper panel) by the final state of model Z)i, an initial condition 
obtained from the parent galaxy and characterized by a diffuse, massive DM halo. The 
positions of this end-product in the {ki,k2) plane are instead remarkably similar to those 
of the one-component models, and the same comments apply. Again, the scatter associated 
with projection effects is smaller than the total thickness of the FP. 

4 ORBITAL ANISOTROPY AND THE FP TILT 

As discussed in the Introduction, a question frequently addressed in the literature is whether 
the so-called FP tilt can be due to some kind of structural and/or dynamical "non-homology" , 
i.e., to a systematic variation from low to high luminosities of the structural and/or dynam- 
ical properties of the galaxies. From an observational point of view, this problem is still in 
general unsettled (see, e.g., Caon, Capaccioh & D'Onofrio 1993, Graham & CoUess 1997, 
Gerhard et al. 2001, Bertin et al. 2001 and references therein) and so here we try to gain 
some hint on its solution by using the results of the numerical simulations described in 
Section 3. In particular we investigate whether relation (7) can be satisfied, in the whole 
observed luminosity range (0.2 < Lb ^40 in the BBF sample), by a systematic variation of 
Ky induced by an appropriate underlying correlation Sa = Ss,{Lb), while maintaining fixed 
the models structure and T*. In practice, the isotropic parent model of each family is placed 
on the FP by selecting its Lb and T^,; then, by increasing its anisotropy and luminosity, 
the family of initial conditions is generated. Note that, at variance with the exploration 
described in Section 3, in this case all the initial conditions are placed by construction on 
the FP, as can be seen from Fig. 5 (upper panel) where initial conditions A2, B2, C2 and D2 
are generated by the parent galaxies A, B, C and D. 

As in the case of Section 3, we found useful to represent the result of the simulations in 
a coordinate system slightly different from the usual k space. In the abscissae axis of Fig. 6 
we plot the quantity k\ — /^g ° that measures how much a given initial condition is displaced 
on the FP from its parent isotropic model, while in the ordinate axis we plot the quantity 
^3 ~ ^3 of the corresponding end-product. As a consequence all the parent isotropic models 
are placed at the origin, and the FP is represented by the solid line k^ = k]^: this line is also 
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Figure 6. End-products of initial conditions placed on the FP at fcg — ^^-j"^ where kij° is the coordinate of their isotropic 
parent galaxy. The dashed lines correspond to Sk^ = a{k3). Note that the abscissae range is significantly smaller than the 
actual range spanned by BBF galaxies, Afca ~ 0.3. Symbols are the same as in Fig. 4. 

the line of stable initial conditions, while the two horizontal dashed lines represent the FP 
thickness. By using an argument similar to that used in the case of Fig. 4, we now obtain 
Sk^ ~ 0.816|A;3 — /cgl and the end-products of unstable initial conditions are vertical segments 
due to projection effects. 

4.1 One— component models 

For each family of 7-models we fixed a constant value of T^, such that the isotropic model, 
with a suitable assigned Lb, lies at the faint end of the FP [ki ~ 2.6, k^ ~ 0.75). Then we 
placed the anisotropic initial conditions on the FP by choosing Lb, according to equation 
(7), as a function of Ky = Kylj, Sg), in order to reproduce the tilt of the FP. 

The inspection of Fig. 6 reveals that it is not possible to reproduce the FP tilt over the 
whole observed range {Aki ^ 2 and Ak^ ~ 0.3 in the BBF sample) by using stable model 
only. In fact, independently of the value of 7, the variations of Ky for Sa^Sas correspond 
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to AA;3<0.04, much smaller than the observed interval. If we consider also unstable (but 
consistent) systems, we can use a wider range of ii'v(7,'Sa) (cfr. Table 1). In this case it is 
possible to reproduce the FP tilt over a much larger interval, Afca ~ 0.2, which is however 
still significantly smaller than the observed one: even if consistency limitations only are taken 
into account, the FP tilt cannot be explained as an effect of a systematic increase of radial 
anisotropy in one-component models, under the assumption of structural homology (see also 
CLR96 and CL97). 

What happens to the end-products of unstable initial conditions is shown in Fig. 6: 
in general they fall well outside the FP thickness, and the departure is larger for larger 
distance from the parent galaxy: they always fall out of cr(/c3) for ki^ — k]^° ~ 0.1: this implies 
that, limiting to end-products of unstable initial conditions which remain inside the FP 
thickness, the maximum amount of the FP tilt that can be explained by pure anisotropy 
is /cg — kf°<0.1. Of course, a larger part of the FP tilt could be covered by considering 
parent galaxies with tangentially anisotropic velocity dispersion tensor (as observed in real, 
low luminosity ellipticals). 

The results presented can also be seen, for a few representative cases, in the upper panel 
of Fig. 5, where the position of the end-products of models A2, B2, C2, D2 is represented by 
the straight line segments. As in the case discussed in Section 3, the end-products remain 
well inside the region populated by real galaxies in the (fci, k2) plane. 

Finally we qualitatively explore an interesting related problem (which however is not 
the argument of our paper), i.e., whether the FP tilt can be reproduced over the whole 
observed range in luminosity, for fixed T^,, in case the requirement of perfect structural 
homology is relaxed. In practice we determined T^, so that the parent isotropic 7 = model 
(to which corresponds the maximum value of Ky) is placed at the faint end of the FP 
(Lb = 0.2, T* = 5). With the same value of T,,, the 7 = 1 and 7 = 2 parent galaxies are 
then placed on the FP, and their positions are characterized by larger Lb (and so ki), due 
to the corresponding decrease of Ky (see first column in Table 1). In this approach, the FP 
tilt can be reproduced over almost all the observed range by using models compatible with 
consistency. However, even in this mixed structural and dynamical non-homology approach, 
if one limits to stable models only, the maximum available range is reduced to Ak^ 0.16, 
approximately half of the required k^ variation. 
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4.2 Two— component models 

In analogy with the one-component case, for each family of (1,1) models, given /i, jS (the 
same as for the models presented in Section 3.2) and T,,,, we placed the parent isotropic 
model at the faint end of the FP and we derived Lb of the initial conditions according to 
equation (7), as a function of Ky = Ky^h, [3, Sg). As already discussed in Section 3.2, for 
each family of (1,1) models we study only two values of the anisotropy radius, s^, = 0.3 and 
Sa = 0.7. We found that for each family these initial conditions, under the only requirement 
of consistency (sa>Sac), span a range Afc3<0.2, smaller than the observed one (see Fig. 6). 
The main difference with respect to the one-component case is that in the {ki,k^) space 
the end-products of the unstable (1,1) models move less than the one-component models 
with similar initial Sk^ and only the most anisotropic systems fall outside the FP thickness. 
An example of these systems is model D2, whose behavior in {ki,k^) and {ki,k2) planes, 
clearly represented in upper and lower panel of Fig. 5 respectively, is similar to that of 
one-component models {A2, B2, 6*2)- Again we found that the DM halo concentration is an 
important quantity, strictly related to the displacement of the end-product with respect to 
the initial conditions. 

These findings suggest that even in the two-component case under the assumption of 
perfect structural homology and constant mass-to-light ratio the FP tilt cannot be explained 
as a consequence of a systematic increase of radial orbital anisotropy. 

5 DISCUSSION AND CONCLUSIONS 

With the aid of numerical simulations of one and two-component galaxy models we explored 
the constraints imposed by the observed thickness and tilt of the FP on the amount and dis- 
tribution of radial orbital anisotropy in elliptical galaxies. The main results are summarized 
below. 

• Remarkably, all the explored models (both one and two-component, and quite indepen- 
dently of the density profile) are found to be unstable when their orbital radial anisotropy 
is high enough to place them outside the observed FP thickness (under the assumption that 
their isotropic parent models lie on the FP). On the contrary all stable models lie inside the 
FP thickness. 

• The end-products of one-component unstable models initially placed outside the FP 
fall back inside the FP: in other words, the larger is the initial displacement from the FP, 

© 2001 RAS, MNRAS 000, 



26 C. Nipoti, P. Londrillo and L. Ciotti 

the stronger is the reassessment of the model structure and dynamics. The behavior of two- 
component models is more varied, due to the fact that the properties of their end-products 
are significantly affected by the amount of mass and the distribution of the DM halo. In 
particular, the end-products of models with massive (either concentrated or diffuse) DM 
halos remain outside the FP thickness, while models with light halos behave essentially like 
one-component models. 

• Since the end-products of the unstable initial conditions are not spherically symmetric, 
their positions on the (fci, ^3) plane depend on their relative orientation with respect to the 
line-of-sight direction. However, the scatter due to projection effects is in general smaller 
than the observed thickness of the FP, both for one and two-component models. 

• We found that it is impossible to reproduce the whole FP tilt with radially anisotropic 
but stable (one and two-component) models under the assumption of constant and 
structural homology. In other words, under these assumptions, luminous galaxies would be 
radially unstable well before the bright end of the FP. 

• At variance with what happens to the end-products of unstable models initially placed 
outside the FP thickness (but exactly for the same reasons), the end-products of unstable 
models with initial conditions on the FP, fall well outside the FP itself. 

Our results lead to some speculation on the formation mechanism and evolutionary sce- 
narios of elliptical galaxies. First, if the (unknown) formation mechanism produces galaxies 
with various degrees of internal radial orbital anisotropy, of which the isotropic ones consti- 
tute the "backbone" of the observed FP, then the most anisotropic systems would be radially 
unstable and would evolve into final states lying on the FP. Then no "ad hoc" fine tuning 
would be required on the amount of radial anisotropy of ellipticals at the moment of their 
formation. In addition, our results concerning the FP tilt could give some indications about 
the importance of dissipationless merging in the history of the assembly of elliptical galaxies. 
In fact , if Es form by hierarchical dissipationless merging, then the very existence of the FP 
necessarily implies structural or dynamical non-homology of the merging end-products. The 
possibility that we explored in this work is that of a substantial dynamical non-homology as 
a function of galaxy luminosity. Our simulations show that this is not a viable possibility to 
reproduce the FP tilt: in this scenario the FP would be destroyed by merging. However, it 
should be clear that we cannot rule out the possibility that merging produces a combination 
of structural and dynamical effects that conspire to maintain galaxies on the FP. For this 
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reason we are now exploring this problem, with the aid of one and two-component galaxy 
merging simulations. 

All our results on the FP thickness and tilt seem to point toward a significant dynam- 
ical homology in real galaxies, and dynamical homology in luminous Es has been recently 
determined by some authors (Gerhard et al. 2001); we note also that an independent ob- 
servational support for dynamical homology is given by the very evidence of the Mbh-o'o 
relation, which relates a dynamical independent quantity (Mbh) with a quantity strongly 
dependent on anisotropy (ctq). The fact that the scatter of the Mbh-c"o is very small means 
that elliptical galaxies are basically dynamically homologous systems. 
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APPENDIX A: INTRINSIC AND OBSERVATIONAL PROPERTIES OF 
THE NUMERICAL MODELS 

In this Appendix we summarize the techniques used to derive the intrinsic and projected 
properties of the end-products of the simulations. 

As discussed in Section 3.1, the intrinsic elhpticities associated with the inertia elhpsoids 
of the end-products of the numerical simulations were used as a measure of their departure 
from spherical symmetry; thus for each numerical simulation we compute the inertia tensor 
lij = Y^k^T'^'^P associated with the density distribution of interest (i.e., stars or halo). 
Following Meza & Zamorano (1997), the sum is extended over all the particles inside r7o, 
the radius of the sphere centered on the center of mass of the galaxy and enclosing the 70% 
of the total mass. The inertia tensor is diagonalized and the ratios of the square root of its 
eigenvalues are used to obtain a first estimate of the elhpticities of the density distribution, 
according to the standard definition 1 — 670/^70 and 1 — 070/070. The orthogonal matrix 
corresponding to the diagonalization is also obtained, and the density distribution is rotated 
accordingly. This procedure is then applied iteratively to the ellipsoid characterized by a = 
r7o and with intermediate and minor axes obtained from the elhpticities computed at each 
stage, up to convergence at some prescribed accuracy level of the elhpticities. Note that 
this procedure automatically select bound particles only: in fact, in all our simulations the 
fraction of "escapers" never exceeds 0.02% of the total number of particles. 

In order to obtain the "observational" properties of our synthetic galaxies, we calculate 
the following projected quantities of their stellar component: the isophotal ellipticity e, the 
circularized effective radius (-R)e, the mean effective surface brightness (/)e, and the central 
velocity dispersion oq. The line-of-sight direction is fixed by the arbitrary choice of d and 
(/?, the two angles of spherical coordinates, expressed in the reference frame where the inertia 
ellipsoid is diagonal. We apply to the system the rotation matrix IZ = lZ2{'d)'Tlz{ip) , where 
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In this coordinate system, the line-of-sight direction coincides with the z axis, while (x, y) 
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is the projection plane. The isophotal axis ratio h/a and the associated elhpticity t = 1 — h/ a 
in the projection plane are determined by using a two-dimensional version of the iterative 
scheme described above. 

The semi- axes and 6e, of the effective isophote are determined under the assumption 
that the ellipticity of the projected density is constant. Finally, the circularized effective 
radius is obtained: 



{R)e = yaebe = aevT^. (A3) 

From the knowledge of (-R)e, the mean effective surface brightness {I) e=L-Q / 2t{ {R) ^ is 
derived modulo the free parameter T,,, = M^./L^. 

The "central" velocity dispersion ctq is computed, by restricting to the central ellipse 
corresponding to a circularized radius {R)e/8. Thus we use 



1 ^ 

where N is the number of particle in the projected ellipse of semi-axes Oe/S and be/8, Vz^ is 
the line-of-sight velocity of the i-th particle, and Vz is the mean velocity integrated along 
the line-of-sight. We find that discreteness effects on the derived values of o"o range from 
0.5% (for 7 = models) to 0.9% (for 7 = 2 models); moreover, the uncertainties on {R)e 
and (/)e never exceed 0.7% and 1.4%, respectively. In other words, the numerical error bars 
on the model measurements are significantly smaller than the observed scatter of the FP 
itself, and so are not important in the context of the paper. 

This paper has been produced using the Royal Astronomical Society/Blackwell Science I^TgX 
style file. 
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